Experimental approximation of the Jones polynomial with one quantum bit 
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We present experimental results approximating the Jones polynomial using 4 qubits in a liquid 
state nuclear magnetic resonance quantum information processor. This is the first experimental 
implementation of a complete problem for the DQC1 model of quantum computation, which uses a 
single qubit accompanied by a register of completely random states. The Jones polynomial is a knot 
invariant that is important not only to knot theory, but also to statistical mechanics and quantum 
field theory. The implemented algorithm is a modification of the algorithm developed by Shor and 
Jordan suitable for implementation in NMR. These experimental results show that for the restricted 
case of knots whose braid representations have four strands and exactly three crossings, identifying 
distinct knots is possible 91% of the time. These results demonstrate the high level of experimental 
control currently available in liquid state NMR. 

PACS numbers: 03.67.Lx, 76.60.-k. 



Quantum information processors have the potential to 
solve some problems exponentially faster than current 
classical methods [T] . While much effort has been concen- 
trated on the most conventional circuit model of compu- 
tation which involves preparation of pure fiducial quan- 
tum states, other models of computation, where only one 
pure quantum bit is required, still offer efficient solutions 
to classically intractable problems. Deterministic quan- 
tum computation with one quantum bit (DQC1) is such 
a model 0. It extracts the power of one bit of quan- 
tum information alongside a register of many qubits in a 
completely random state. Study of DQC1 was originally 
motivated by liquid state nuclear magnetic resonance 
(NMR), which is a high temperature ensemble model 
of quantum computation. Although this model of com- 
putation is weaker than conventional models with many 
pure qubits, it has been shown to have several impor- 
tant applications where classical methods are inefficient: 
simulating quantum systems [2], estimating the average 
fidelity decay under quantum maps [3 , and quadratically 
signed weight enumerators [4] . Additionally, the approx- 
imation of the Jones polynomial at the fifth root of unity 
has recently been shown to completely encapsulate the 
power of DQC1 |5j. DQC1 algorithms have been exper- 
imentally implemented in optics [6 and liquid and solid 
state NMR [H |3] - none of which has been shown to be 
DQCl-complete. In [8], the authors implement a DQC1 
algorithm on two qubits to evaluate the Jones polynomial 
at various points for specific knots. This letter describes 
the implementation of an instance of a DQCl-complete 
algorithm [5], which scales for any size knot. 

Unlike its name suggests, DQC1 does not require a 
completely pure qubit to provide an advantage over 
known classical methods, but rather a small fraction of 
a pure qubit. This pseudo-pure state is almost com- 
pletely mixed with a small bias towards the ground state, 



and is used as the control qubit in the DQC1 algorithm. 
A unitary is performed on the qubits in the completely 
mixed state and is controlled by the pseudo-pure qubit. 
Measurements of (a x ) and {a y ) yield the real and imagi- 
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FIG. 1: The DQC1 circuit where the pure qubit has a bias 
of e towards the ground state. Measurements of (a x ) and 
(<7 y ) will yield the real and imaginary parts of £Tr([/„)/2", 
respectively. 

nary parts of the trace of the unitary, normalized by the 
amount of polarization on the pure qubit. 

Applications for the Jones polynomial are extensive in 
physics; for example, the fields of statistical mechanics, 
quantum field theory, and quantum gravity would benefit 
from an efficient method for approximating this polyno- 
mial [TO]. Knot invariants help to solve a fundamental 
problem in knot theory: determining if two knots, defined 
as the embedding of the circle in M 3 , are topologically 
different, up to ambient isotopy. Two knots can only be 
confirmed identical if one can be maneuvered into the 
other by a sequence of Reidemeister moves, which keep 
the topological properties of knots intact. This process is 
very tedious as often the sequence of Reidemeister moves 
require an increase in the number of crossings in the knot. 
Even the simplest such problem of identifying the unknot, 
a circle with no crossings, has been shown to be contained 
in the complexity class NP [TT]. Knot invariants, such as 
the Jones polynomial, have the same value for different 
representations of the same knot. In other words, if a 
knot invariant evaluates to different values for two knots, 
they are guaranteed to be distinct. This makes them a 
welcome alternative to sequences of Reidemeister moves. 
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Unfortunately, exact evaluation of the Jones polynomial 
at all but a few points is #P-hard [12]. Several efforts 
for finding quantum algorithms for the Jones polynomial 
have been attempted and approximations at several spe- 
cial points have been shown to be BQP-complete [T31I14) . 
Largely building on this work it was then shown that ap- 
proximations of the Jones polynomial for trace or plat 
closures at principal roots of unity can be computed on 
a quantum computer in polynomial time [15 . Later it 
was shown that for the plat closure the problem is BQP- 
complete [TSJ [T7] . The algorithm developed by Shor and 
Jordan shows that approximating the Jones polynomial 
at the fifth root of unity for any knot is a complete prob- 
lem for DQC1. 

For the purposes of this algorithm, knots are described 
in the discrete language of braid groups. Every knot can 
be written as a braid, which is a series of strands crossings 
over and under each other with loose ends at both the 
top and bottom. Braids can then be converted into a 
knot by the trace closure which connects the top and 
bottom ends of the braid in sequential order. The braid 
group for m strands B m is generated by s\ . . . s m _i that 
denote elementary crossings where Si indicates the i th 
strand crossing over the (i + l) th strand and s^ 1 indicates 
strand (i + 1) crossing over strand i. These elementary 
crossings satisfy the relations: SiSj = SjSi for \j — i\ > 1 
and s l+ is. ( s i+ i = s i s i+1 s l . 

The implemented algorithm utilizes the Fibonacci rep- 
resentation of the braid group B n , which is described in 
the context of the Temperley-Lieb recoupling theory [IB] . 
In this theory there are two particles p and *, which ex- 
hibit the following properties: p interacts with another p 
to create apora* particle, * interacts with a p to always 
create a p particle, and two *'s never interact. Strings of 
these particles create a basis in a complex vector space. 
More details of this representation can be found in [TS] . 
but for our purposes it suffices to state that for a braid 
with m strands the basis vectors contain m+1 elements of 
p's and *'s with the restriction that no two * particles be 
beside one another. These basis vectors are then trans- 
formed into the computational basis and unitary matrices 
<7j, which represent each elementary crossing in the braid 
group, are constructed. For the particular form of these 
unitaries, please refer to [5]. 

The algorithm developed by Shor and Jordan approxi- 
mates the Jones polynomial at the single point t = e 2 "/ 5 
by finding the weighted trace of a unitary that describes 
the braid representation of the knot. The algorithm is 
modified for this implementation and the varied por- 
tions are described below. The primary difference is 
in the encoding of the basis states. The Fibonacci ba- 
sis vectors consist of four distinct subspaces, only two 
of which arc relevant for the algorithm: the /,„ vec- 
tors of the form and / m _i of the form 
where /„ = [1, 1, 2, 3, . . .] is the fibonacci sequence. 
These are the only two subspaces that are encoded in 



this implementation. The Zeckendorf representation, 
z' = 2"~ 1 si + J2i^2 s i+sifi converts the Fibonacci ba- 
sis vectors into integers that are then converted to a 
non-saturated computational basis. The second notable 
difference is the method used to calculate the weighted 
trace, defined as 

WTr = 1 x (trace of subspacc *...*) + 
<f> x (trace of subspace * . . . p) , 

where (f> = (1 + v5)/2 is the golden ratio. Implement- 
ing these weights for our encoding is achieved by purify- 
ing the second qubit, then applying a rotation taking |0) 
to (-^010) + |l))/-^/l + (j>, which ensures that each sub- 
space receives the desired weight. The computational 
model now contains two initialized qubits, however this 
modification does not change the computational power 
as DQC(fc) is known to have the same computational 
power as DQC1 for k that can grow logarithmically with 
the total number of qubits [5]. The extra basis states 
are accounted for in the final calculation of the Jones 
polynomial. The circuit for our evaluation of the Jones 
polynomial for braids with four strands can be seen in 
figure [2j It is worthwhile to note that the off-diagonal cl- 
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FIG. 2: Circuit diagram for the approximation of the Jones 
polynomial for the knots whose braid representations consist 
of four strands. The initial state given is the traceless devia- 
tion matrix. The single qubit gates are the Hadamard and the 
rotation for implementing the weights of the trace. The mea- 
surements performed on the top qubit are expectation values 
of the Pauli-x and y operators. 

cmcnts in the rotated pure qubit do not contribute to the 
algorithm as the unitary matrices U n are always block di- 
agonal, thereby eliminating the off-diagonal elements in 
the calculation of the trace. The state of the top qubit 
at the completion of the algorithm is 

1 / 1 WTr([/t) \ 

1 2»- 1 (l + ^)U Tr ( [/ ») 1 /' 

which upon measurement of (a x ) and (a y ) yields the real 
and imaginary parts of M. = WTr(£/ rl )/(2 ,l ~ 1 (l + </>)) re- 
spectively, where n is the number of qubits in the bottom 
register. The measured quantity M. is then used to cal- 
culate the approximation of the Jones polynomial, V(t), 
corresponding to the trace closure of the given braid at 
t = e i2 "/ r °, 

y( e <2*/6) = (_( e i2,r/5)4)3^-l (V-^l + (j>)M - Kj , 
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where re = (2™" 1 — f m )4> + (2" _1 — f m -i) and w is the 
writhe of the braid, defined as the number of positive 
crossings minus the number of negative crossings. 

Liquid state NMR offers one of the most advanced im- 
plementations of quantum information processors with 
high fidelity control of multiple qubits [T^]. The qubits 
are a bulk ensemble of identical spin- 1/2 nuclei that ex- 
hibit a two-level energy structure in the presence of a 
strong magnetic field. The ensemble of approximately 
10 20 molecules are manipulated in parallel and an en- 
semble measurement is performed using quadrature de- 
tection of the free induction decay to give (a x ) and {<r y ). 
The algorithm described above was demonstrated in liq- 
uid state NMR for the set of knots whose braid repre- 
sentations have four strands and three crossings. There 
are six distinct knots in this set and hence, six distinct 
Jones polynomials. The goal of the experiment is to dis- 
tinguish between two distinct knots given their braid rep- 
resentations. The subspaces of interest have fn = 3 and 
fi-i = 2 basis states respectively, thus the encoding of 
the basis states requires 3 qubits in the bottom register 
and a fourth as the control qubit. 
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FIG. 3: The molecule trans-crotonic acid and a table with the 
parameters of the Hamiltonian given in Hz. The shaded diag- 
onal elements represent the chemical shifts Ui with the Hamil- 
tonian Ei7rcjiCT*. The remaining elements indicate the scalar 
coupling constants Jij with the Hamiltonian 'Ei < jnJijala 3 z . 

The experiment was implemented on a Bruker Avance 
700 MHz spectrometer using the molecule trans-crotonic 
acid (shown in fig [3| . The four qubits are experimen- 
tally realized by the four carbon nuclei, synthesized to 
be carbon-13, while the hydrogen are decoupled using 
the WALTZ-16 [2U] composite pulse sequence. C\ is 
our readout qubit whose initial state is the thermal 
state of p = 1 + eZ, Ci is purified to the pseudo-pure 
state |0)(0|, and the remaining C3 and C4 are initial- 
ized to the completely mixed state. The radio frequency 
(r.f.) pulses that implement the unitary transforma- 
tions are numerically generated using the GRAPE al- 
gorithm [211 [22] which starts from a random guess and 
is then iteratively improved through a gradient ascent 
search. The GRAPE pulses are optimized to produce a 
fidelity \ti'(Ug 0al U S i m )\ 2 / d 2 , where d is the dimension of 
the Hilbert space of U goa i, of no less than 0.998 and are 
designed to be robust to small inhomogeneities (±3%) in 
the r.f. control field. Each controlled-cjj unitary transfor- 
mation is designed as a single pulse of 60 ms. The pulses 
are corrected for non-linearities in the pulse generation 



and transmission to the sample by measuring the r.f. sig- 
nal at the position of the sample using a feedback loop 
and iteratively modifying the pulse accordingly. Through 
the feedback loop the implemented pulse can be mea- 
sured and was found to have a simulated fidelity of 0.99 
after correction. 

The resulting spectrum is fit and compared to a refer- 
ence spectrum, traditionally of the initial state, to give 
the expectation value results. In this experiment, pulses 
whose propagator was designed to be the identity were 
generated using GRAPE to have the same length and 
the same average power and fidelity as the controlled- 
Oi. These pulses were implemented and used to create 
a reference spectrum in an attempt to normalize some 
decoherence effects. The state measured after three suc- 
cessive identity pulses, totaling 180 ms had only 60% of 
the original signal (see figure [I| , indicating this as a cru- 
cial step in the experimental procedure. 
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FIG. 4: (Color Online) The dashed (red) spectrum represents 
the pseudo pure state immediately after creation in the top 
(a) graph. The solid (black) spectrum is the same pseudo 
pure state after 180 ms of pulses designed to perform the 
identity. In the bottom graph (b), the solid spectrum (black) 
indicates the final state of the experiment and it is compared 
to simulation (dashed (blue)). This particular experiment is 
for the knot whose braid representation has crossings S1S2S3. 

The algorithm was implemented for 18 different braids, 
which correspond to 6 distinct Jones polynomials. The 
results are displayed in figure [5j Systematic errors from 
imperfect initial state preparation and decoherence not 
captured by the reference state result in the offsets from 
the theoretical values. The main contribution to the 
spreading of the experimental points is the finite fidelity 
of the optimal control pulses. 
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Approximation of the Jones Polynomial 
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FIG. 5: The results for the approximation of the Jones poly- 
nomial for knots whose braid representations have 4 strands 
and 3 crossings. There are six unique knots of this kind and 
their theoretical values of the Jones polynomial are plotted 
for the six experiments. The corresponding experimental data 
points of three braid representations for each experiment are 
plotted along with error ellipses demonstrating the statistical 
error (with 86.5% confidence levels or 2a). The distribution is 
generated by simulating each experiment 200 times with sin- 
gle pulse fidelities of 0.99 which is the implemented pulse fi- 
delity. Using the error ellipses as discriminators, these results 
yields a 91% success rate for distinguishing distinct knots. 



Two values of the Jones polynomial at best can distin- 
guish between two knots if they are sufficiently far apart, 
and at worst, give no information, as even evaluations 
of the Jones polynomial that are identical would not be 
sufficient information to conclude the two knots are iden- 
tical. This leads to two types of errors when interpreting 
the data: passive and fatal errors. Passive errors occur 
when two distinct knots are impossible to distinguish be- 
cause of their relatively close distance to one another, 
while fatal errors occur when two identical knots are de- 
termined to be distinct. The success rate for determining 
whether knots are distinct is calculated as the average 
of the percent of distinct knots correctly identified and 
the percent of identical knots correctly indistinguishable. 
The error ellipses give a direct method for determining if 
two knots are distinct. If the error ellipses for a pair of 
knots do not overlap then it is inferred that the knots are 
distinct, whereas if the two ellipses overlap no informa- 
tion is gained. For the confidence region plotted in figure 
[5} 134 of the possible 135 pairs of distinct knots are cor- 



rectly distinguished with 3 fatal errors of a possible 18, 
corresponding to a success rate of 91%. 

Approximation of the Jones polynomial is an example 
of a classical problem that appears intractable, but that 
can be solved using one clean qubit quantum computers. 
This is the first experimental implementation of a DQC1- 
complete problem, and is performed in liquid state NMR 
with four qubits, resulting in a 91% success rate for braids 
with four strands and a total of three crossings. In future 
work it will be interesting to see how the values of the 
Jones polynomial spread as you scale to larger knots and 
what size knot can be experimentally implemented before 
noise and control errors destroy the quantum advantage. 

G.P. would like to thank M. Ditty for his technical 
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by NSERC, Quantum Works, and CIFAR. 
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